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Abstract 

We have analyzed data on the free ion yield observed in liquid isooctane irradiated by 
60 Co 7 rays within the framework of the Onsager theory about initial recombination. 
Several distribution functions describing the electron thermalization distance have 
been used and compared with the experimental results: a delta function, a Gaussian 
type function and an exponential function. 

A linear dependence between free ion yield and external electric field has been 
found at low electric field values {E < 1.2 MV/m) in excellent agreement with the 
Onsager theory. At higher electric field values we obtain a solution in power series 
of the external field using the Onsager theory. 
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1 Introduction 



Liquid filled ionization chambers are currently used in radiotherapy both for 
dosimetry ([1], [2], [3], [4]) and portal imaging [5]. One of the most commonly 
used liquids is isooctane (2,2,4 trimethylpentane). This nonpolar liquid has a 
quite constant stopping power ratio to water in a very wide energy spectrum 
(less than 3% variation from 0.1 MeV to 20 MeV [1]) and also its intrinsic 
mass density allows to achieve a spatial resolution in the millimeter range for 
therapy beams [5]. 

Free ion yield Gfi(E,T), is defined as the number of electron-ion pairs 
escaping initial recombination per 100 eV of absorbed energy, and experi- 
mentally it is obtained from ionization current measurement. The knowledge 
of how it varies with temperature T, with external electric field E, and with ra- 
diation type, constitutes a fundamental problem to understand the operation 
of these devices. These dependences have been measured in a large number 
of liquids ([6], [7]). The Onsager theory [8] describes Gfi(E,T), and has been 
tested in several liquids with good results (see for example [7]). The Onsager 
theory predictes a linear relationship between ionization current and electric 
field at low electric field values. The previous dependence can be obtained 
from numerical resolution of the Onsager theory. This linear behavior has to 
be extrapolated to very low electric field strength region because volume re- 
combination depletes free charge density produced by radiation in the liquid. 

In the current work we describe a detailed method to apply and to test the 
Onsager theory at the low field strength region and to obtain a precise depen- 
dence of Gfi with electric field in the linear region for liquid isooctane. 



2 Theoretical considerations 

2. 1 Onsager theory of initial recombination 

When ionizing radiation interacts with a liquid, electrons released from 
molecules thermalize at a distance r, where electron and positive ion are still 
bounded by the Coulomb interaction. This will cause the recombination of 
the primary ionization pairs produced, which is called initial recombination. 
These effects are much more relevant in liquids than in gases due to the fact 
that mass density of liquid hydrocarbons is almost three orders of magnitude 
higher than density of gases at normal conditions. 

Onsager solved the problem of the Brownian movement of an electron under 
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the influence of both the ion Coulomb attraction and an external electric field 
E [8]. For isolated ionizations, initial recombination escape probability of an 
electron-ion pair within the Onsager theory is 

Er 

$(r, E, 0, T) = exp{- — — (1 + cos 6)} 

EqT o 
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r /r 

x exp(-s) ds (1) 



where r is the initial separation between electron and ion (i.e. the thermal- 
ization distance), G is the angle between the line that initially connects the 
electron-ion pair and the external electric field. The variables r = e 2 /4ireK,T 
and E = 2/tT/er are the Onsager radius (the distance at which Coulomb 
energy equals thermal energy kT) and the Onsager field (the field that would 
produce a voltage 2nT/e over a distance r ). Here e is the liquid dielectric 
constant (e = 1.94 • eo for liquid isooctane at room temperature), T is its tem- 
perature and k is the Boltzmann constant. Finally, Jo denotes the zeroth-order 
Bessel function. 

Mozumder [9] converted the integral of equation (1) in an infinite series using 
properties of the Bessel functions. He also eliminated the angular dependence 
averaging over a uniform distribution of cos 0. Then, the angle averaged escape 
probability takes the following form, 

Tp r CO r> rp 

%V) = 1-^E^ )Mr«M (2) 



where A n (x) is the n order incomplete gamma function, which is given by: 

oo ^.k—n oo ^.m 

A n (x) = exp(-x) Yl (k „n| = exp(-x) ]T — 

fc=2n+l \ K ~ n ) 1 m=n+l m! 

= 1 " exp(-x) [ 1 + x + X - + • • • + X — ] (3) 



The next expression is more practical for numerical computation of U(r, E, T): 

A n+1 (x) - A n (x) = -(x n+1 /(n + 1)!) exp(-x) (4) 
Aq(x) = 1 — exp(-x) (5) 
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Equation (2) is the most adequate formula for calculating escape probabilities 
for arbitrary values of initial separation and external electric field. In fact, 
it will be the formula that we will use for numerical calculations of escape 
probabilities. The expansion in power series of the external field is implicit in 
equation 2. If we operate in equation (2) then we can obtain it explicitly: 



U(r,E,T) =exp(-r /r) 



oo , rp 
n=l V " t/ 



(6) 



where B n (x) is a polynomial of order n — 1 in x, which takes the following 
form: 



Bn{x) = £ 



m=l L \k 



x 
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The numerical coefficients Fj? were calculated by Mozumder [9] , and they are 
given by 



and 



F£ =0 for k>n 
2 n 



(n + l)\ 
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for k < n 



(10) 



We must keep in mind that thermalization distance is not the same for all elec- 
trons. Due to this fact a distribution function /(r), such as to J °° f(r)dr = 1, 
is usually introduced ad hoc to describe electron thermalization distances. Sev- 
eral distribution functions had been used by different authors in the literature 
(see for example [7] and [9]). In the current article we will test the three most 
used: 

• delta function 



f(r,p) = 5(r- p) 
Gaussian type function 



(11) 



f(r,P,cr) 



iVexp[-(r - p) 2 /(T 2 ] ifr>0 
if r < 
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exponential function 

- exp(— r/p) if r > 



fir, p) 



(13) 

if r < 



The first and second distributions are characterized by a single parameter p. 
The Gaussian type distribution, in addition, has a dispersion parameter a. 
However, in order to obtain a single parameter distribution we will correlate 
both parameters. In advance we will take a = 0.25 • p such as [9]. This choice 
provokes that the difference between the normalization factor N, and the 
normal Gaussian normalization factor l/x/^er, be less than 0.1%. 

Within this framework we can write the escape probability averaged over the 
thermalization distances as: 

oo 

P esc (E, T) = J n(r, E, T) f(r) dr (14) 



An interesting property of the Onsager series is that the first expansion term 
of the escape probability does not depends on the thermalization distance 
r. From equations (6) and (7) we can derive that when the electric field is 
sufficiently low to verify 

3(§)»|l-2(-)| 
E r 

the n — 1 term is much higher than the n = 2 term in equation (6) (and hence 
than terms n > 2). Then we can truncate the series to first order and so the 
escape probability rises linearly with the electric field: 

U(E, T) = exp(-r /r) (l + + • • •) (15) 
In this case the intercept-to-slope ratio E c is predicted to be the Onsager field 

Eq. 

Ec = Eo = s_H^f 



This result is a powerful test to compare experimental data with the Onsager 
theory because it does not depend on the distribution function used to describe 
the electron thermalization distance. 
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Fig. 1. Escape probability calculated using a Gaussian distribution (continuous line), 
an exponential distribution (dotted line) and a delta distribution (dashed line), plot- 
ted against electric field for liquid isooctane at T = 294 K. 
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Fig. 2. . Escape probability calculated using a Gaussian distribution (continuous 
line), an exponential distribution (dotted line) and a delta distribution (dashed line), 
plotted against electric field for liquid isooctane at T = 294 K. We can see how the 
relationship becomes linear at low electric fields (approximately E < 1.2 MV/m for 
Gaussian type and delta distributions, E < 0.5 MV/m for exponential distribution). 

As we will see later, for isooctane at room temperature irradiated by 7 photons, 
r/r ~ 0.6 and the linear approximation (15) is right at low electric fields. 
Figure 1 shows escape probability variation with external electric field for the 
case of liquid isooctane at T = 294 K. Figure 2 shows the low electric field 
region where a linear relationship is expected. 
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2.2 Free ion yield calculation within the Onsager framework 

Free ion yield Gfi, is defined as the number of electron-ion pairs created (i.e. 
those escaping initial recombination) in the ionization medium per 100 eV 
of absorbed energy. This magnitude plays a similar role to the W factor in 
gases. However, W is constant and free ion yield depends on temperature, on 
external electric field and on radiation type as escape probability does. Within 
the Onsager theory we can write 

G fi (E,T)=N tot P esc (E,T) (17) 

where N tot is the total number of electron-ion pairs formed initially in the 
ionization medium (i.e. before initial recombination) per 100 eV of absorbed 
energy, and P esc is the escape probability described in the previous subsection. 

Free ion yield at zero external electric field is denoted as G ^, and takes the 
following form: 

oo 

G fi = N tot J f(r) exp(-r /r) (18) 
o 

If we are working at electric field values for which the relation between escape 
probability and electric field is linear, then we can write 

= G° fi + aE (19) 

and the intercept-to-slope ratio from equation (16) 

E C = ^ = E (20) 



3 Experimental results 

In order to measure the ionization current from an isooctane layer under irra- 
diation we built a square shaped parallel plate liquid ionization chamber. The 
chamber walls were fabricated using FR4 fiber glass reinforced epoxy copper 
clad on both sides, covering a total area of 4.08 cmx4.72 cm. The FR4 thick- 
ness was 0.8 mm while the copper layer was 35 /im thick. The two chamber 
walls were glued on both sides of an epoxy plate spacer to provide the 0.8 mm 
isooctane gap. To guarantee a constant dose rate along the gap, the detector 
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Fig. 3. Scheme of the liquid ionization chamber cross section. 

was inserted between two PMMA 2 plates of 5 mm and 5.8 mm thickness, 
the first one on the incident beam side. As ionization medium we used Merk's 
liquid isooctane 3 with an estimated purity of 99.8%. Figure 3 shows a scheme 
of the device cross section. 

Experimental tests of the chamber were made in the 60 Co unit of the Complexo 
Hospitalario Universitario de Santiago (CHUS). The radiation field was set 
to 10 cmxlO cm at the isocenter 4 to cover the whole chamber. In order to 
obtain isooctane free ion yield we measured the ionization current produced in 
the whole chamber using a nanoammeter Phillips Fluke PM 2525, for several 
polarization voltages. High voltage was supplied by a CAEN N471A NIM 
module. 

Distance between the detector and the cobalt source was set to 130 cm (equiv- 
alent to a dose rate around 0.4 Gy/min). The dose rate was chosen to have 
a negligible volume recombination in the upper part of the ionization current 
vs. voltage curve. Figure 4 shows the experimental data obtained. 

If volume recombination and space charge effects can be ignored, the ionization 
current is proportional to the number of electron-ion pairs released in the liquid 
per unit time and unit volume N in (initial recombination is included in N in ): 



I = ehAN ir 



(21) 



2 Polimethylmetacrylate. 

3 Isooctane Merk Uvasol quality grade. 

4 80 cm from the cobalt source. 
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where e is the electron charge, h is the isooctane gap and A is the detector 
area. From this equation the free ion yield can be calculated as 



G f , = ^ (22) 



where Ae is the energy deposited in the medium per second. In this case 
Ae = (4.79 ± 0.19) • 10 11 100 eV/s, that was calculated through the EGSnrc 
code. 

To apply equation (22) we require a charge collection efficiency higher than 
99% and also that field screening be negligible. To calculate this limit we used 
a numerical simulation of the charge carriers transport. When the distance 
between the cobalt unit and the detector is 130 cm, and the polarization 
voltage is higher than 600 V, general charge collection efficiency is higher 
than 99%. This agrees with the Greening theory ([10], [11]) about general 
charge collection efficiency. Within this theory the polarization voltage that 
must be applied to obtain an efficiency /, is 

V 2 = n 23 

6 Cj-1) 1 ; 



with 



a 



m 



ck-^. k_ 



In these equations k + and k_ are the mobilities of the positive and negative 
charge carriers and a is the volume recombination constant (we used k + = 
k_ = 3.2 • 10 -8 m 2 V _1 s _1 and a — 5.9 • 10~ 16 m 3 s -1 ) 5 . Introducing the 
numerical data of the experimental set-up in equation 23 we obtain / > 0.99 
for V >625V. Then we only apply equation (22) to data for V >600 V. 

For lower voltages the free ion yield has to be extrapolated because volume 
recombination and space charge effects provoke charge losses and field screen- 
ing, and equations (21) and (22) no longer hold. We expect from section 2.1 
a linear relationship between the free ion yield and the electric field at low 
electric field values. Figure 4 shows this linear relationship between 600 V and 
1000 V. At higher voltages the relationship begins to deviate from linearity 
as shows the figure. Then, we use data in the range 600 V < V <1000 V to 
extrapolate the free ion yield at low electric field region. 



5 Values obtained under irradiation with X rays, in agreement with [12]. 
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Fig. 4. Ionization current measured in the chamber against polarization voltage. The 
continuous line shows the linear extrapolation at low field strength region. 

This experimental linear relationship between free ion yield and external elec- 
tric field E, is 



Gfi(E) = (0.32 ± 0.02) + (1.73 ± 0.09) • 10" 7 • E 



(24) 



where E is given in V/m and the free ion yield in pairs/100 eV. Equation (24) 
is valid up to E = 1.2 MV/m with a confidence level of 96%. For higher values 
be have to take into account more terms in the equation (6). 

For isooctane the total number of electron-ion pairs produced per 100 eV 
of absorbed energy is N tot = 1.83 (also calculated with the EGSnrc code). 
Inserting this value and the obtained G ^ in equation (18 ) we can obtain 
the parameter p, for the delta (11), Gaussian type (12) and exponential (13) 
distribution functions. The results are p = 168 ± 6 A for first and second 
distributions, and p = 178 ± 10 A for the exponential one. 

Taking the numeric values for N to t and p we can calculate the theoretical pre- 
diction for free ion yield within the Onsager theory using equations (17), (14) 
and (2), and compare these theoretical results with experimental data. Fig- 
ure 5 shows results obtained using the three considered distribution functions. 
Delta and Gaussian type distributions agree with the experimental data, but 
not the exponential distribution. 

The intercept-to-slope ratio from equation (24) is E c = (1.83 ± 0.12) MV/m. 
Within the Onsager framework we obtain (see equation (16) and figure 2), 
E c = E = (1.74 ± 0.02) MV/m at the room temperature, T = (294 ± 2) K. 
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Fig. 5. Free ion yield against external electric field in the linear region. The diamond 
points show the experimental points for E>0.75 MV/m (for lower fields we have to 
extrapolate using equation (24))- The dotted, dashed and continuous lines correspond 
to the theoretical prediction using exponential, delta and Gaussian type distribution 
functions respectively. 

Experimental value and theoretical value for E c are clearly in agreement. 



4 Conclusions 

We have obtained and analyzed data on isooctane free ion yield irradiated 
by 7 photons, from a cobalt source, within the framework of the Onsager 
theory. Three distribution functions (describing separation distance between 
electron-ion pairs when thermalization is achieved) have been considered: a 
delta function, a Gaussian type function and an exponential function. The 
first and the second describe data correctly in the covered electric field range, 
but not the exponential function. 

This fact means that free ion yield depends in a fundamental way on the choice 
of the distribution function, which is not predicted by the theory. The good 
agreement between the experimental data and the theoretical prediction using 
a delta or a Gaussian type distribution with a dispersion parameter a = 0.25 -p 
seems to show that in this case /(r) is a Gaussian type function with a small 
dispersion parameter. 

If electron would suffer a large number of independent collisions before ther- 
malization, a Gaussian distribution function for the thermalization distances 
will agree with the central limit theorem. The lack of data about the nonpolar 
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liquids cross sections makes difficult to obtain models describing the nature 
of f(r). Some computer simulations have been made (see for example [13]) 
in this sense, however more theoretical and numerical work is needed in this 
area. 

On the other hand, the theoretical prediction for the intercept-to-slope ratio 
E c , is in agreement with the experimental value. This is an important result 
to check the Onsager theory because it does not depend on f(r). 
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